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RESEARCH  SUMMARY 


THE  AUTHOR 


The  original  assumption  (made  by  Rothermel  in  his 
1972  paper— A  Mathematical  Model  for  Predicting  Fire 
Spread  in  Wildland  Fuels)  that  the  reaction  intensity  in 
the  combustion  zone  is  linear  with  loading  is  experi- 
mentally supported  within  the  fuel  sizes,  loadings,  and 
packing  ratios  commonly  encountered  in  wildland  fuels; 
further  confirmation  is  needed  in  light  loadings  of  the 
fine  fuels. 

Rothermel's  empirical  relations  for  the  propagating 
flux  ratio  are  judged  to  be  adequate  for  their  current  use, 
though  perhaps  they  underestimate  the  propagating  flux 
in  heavy  loads  of  fine  fuels.  New  evidence  is  presented  to 
physically  relate  the  propagating  flux  ratio  to  the  optical 
density  (a/3)  of  the  fuel  bed.  Again,  further  development  is 
needed. 

Experimental  difficulties  in  rigorous  determination  of 
reaction  zone  size  remain  unsolved.  Thus,  physical  inter- 
pretation of  reaction  intensities  should  be  made  with 
caution.  Similarly,  the  causes  and  consequences  of  the 
optimum  packing  ratio  should  be  interpreted  with 
caution;  for  example  (in  the  author's  opinion),  the  idea 
that  reaction  intensity  is  limited  by  oxygen  supply  (in 
draft)  at  heavier  packing  ratios  is  premature  and  has 
been  overemphasized. 
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INTRODUCTION 

A  mathematical  model  for  fire  in  wildland  fuels 
(Rothermel  1972)  is  the  basis  of  several  forest  fire 
management  systems.  During  development  of  the  model, 
Rothermel  found  it  necessary  to  make  several  assump- 
tions about  the  effects  of  fuel  geometry  and  fuel 
moisture  upon  combustion.  The  experimental  work  de- 
scribed here  was  initiated  to  examine  these  assumptions 
and  to  extend  the  data  base  beyond  the  range  of 
Rothermel's  early  work.  The  interactions  between  fuel 
moisture,  load,  depth,  packing  ratio,  and  fuel  particle  size 
are  of  primary  concern,  because  some  combinations  of 
these  variables  represent  fuel  arrays  that  will  not  support 
fire.  We  consider  only  the  no-wind  and  zero-slope  case. 

Analysis  of  data  from  more  than  250  experimental  fires 
indicates  clearly  that  additional  data  (probably  100  fires) 
are  needed  before  definitive  corrections  can  be  made  to 
the  fire  model  equations.  Although  the  results  presented 
here  are  preliminary  and  tentative,  some  of  the  original 
assumptions  are  confirmed:  the  original  functional  form 
of  the  spread  equations  is  retained  for  the  geometric 
variables  (o.  p.  6.  Wq.  etc.)  while  the  numerical  coefficients 
are  reevaluated;  a  simplified  fuel  moisture-damping 
coefficient  is  proposed:  a  separate  probability  function  is 
introduced  to  predict  extinguishment.  Reaction  and 
critique  are  solicited. 


PRELIMINARY  DISCUSSION 

Following  the  lead  of  Frandsen  (1971),  Rothermel  (1972) 
developed  the  basic  spread  rate  equations  based  on 
physical  principles.  One  can  begin  with  the  proposition 
that  the  propagating  flux  (Ip)  required  to  drive  a  fire  is 
proportional  to  the  heat  needed  to  bring  the  fuel  to 
ignition  (p^£Q.g);' 

Ip=R(Pb^Qig)  (1) 

where  the  constant  of  proportionality  is  the  spread  rate, 
R.  The  preignition  term  (  p^^Q  )  is  experimentally  deter- 
mined by  calorimetric  methocTs  and  is  not  considered 


here.  The  propagating  flux  is  dependent  on  the  fuel 
particle  size,  arrangement,  chemistry,  and  environmental 
variables. 

We  next  postulate  that  the  propagating  flux  is  propor- 
tional to  the  reaction  intensity,  Ip  (the  power  density 
released  by  the  fire), 

Ip  =  4Ir  (2) 
and  assume  that  I  is  a  function  of  only  fuel  bed  geometry 
(i.e.,  the  fraction  of  energy  transferred  to  the  fuel  is  a 
function  of  particle  size  and  packing  ratio,  etc.). 

An  important  key  to  Rothermel's  success  is  the  con- 
cept of  independently  evaluating  Ip  and  Ip  so  that 
spread  rate  can  be  functionally  related  to  the  reaction 
intensity  defined  by: 

Ip  =  -h(dw/dt).  (3) 

Then  by  a  series  of  definitions  and  logical  assumptions 
the  reaction  intensity  (Rothermel's  equation  27)  becomes 
a  product  of  mathematically  separable  functions  of  the 
primary  independent  variables  (o,  p,  6,  w^,  M,,  etc.) 

Ir  =      hr>^  (4) 
In  summary,  it  has  been  assumed  (heretofore)  that 

•  The  propagation  flux,  Ip,  is  linearly  proportional  to 

reaction  intensity,  Ip; 

•  The  proportionality  constant,  |,  depends  only  on 

fuel  particle  size  and  fuel  bed  packing  ratio; 

•  Ip  is  linearly  proportional  to  net  loading,  w^; 

•  The  moisture-damping  coefficient,  r^^,  is  a 

function  of  fuel  moisture,  m,,  alone; 

•  The  mineral-damping  coefficient,  ^s-  depends  only 

on  the  silica-free  mineral  content  of  the  fuel: 
and 

•  Optimum  reaction  velocity  (the  average  net  frac- 

tional dry  weight  loss  rate  in  the  combustion 
zone)  is  a  function  of  packing  ratio  only  for 
fuel  beds  with  only  one  size  of  fuel  particle,  a. 


'  The  nomenclature  is  defined  in  table  1.  "Basic  Fire  Spread 
Equations."  reproduced  in  metric  form  from  p.  26.  27.  Rothermel's  paper, 
and  Wilson  (1980). 
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Furthermore,  the  empirical  damping  coefficients, 
and  rjg,  are  defined  such  that 

•  ^7i^  =  1  at      =  0  and 

•  f7g  =  1  at  Sg  =  0;  thus, 

•  Ip  =      h  r'  properly  describes  the  fire  state 

in  the  limit  of  no  fuel  moisture  and  no 
minerals. 

As  shown  below,  the  experimentally  determined  reac- 
tion intensity,  Ip,  is  dependent  on  the  experimentally 
determined  reaction  zone  size,  and  the  propagating  flux 
ratio,  I,  is  determined  as  the  fraction  of  Ip  that  is 
required  to  drive  the  fire  at  the  measured  spread  rate,  R. 
Thus,  any  consistent  experimental  technique  for  deter- 
mining the  reaction  intensity  will  provide  a  consistent  set 
of  equations  for  Ip  and  |  that  will  predict  the  spread  rate, 
R  (within  the  range  of  experimental  applicability).  The 
reaction  intensity,  however,  is  also  used  in  applied 
models  to  predict  other  fire  phenomena  such  as  total 
heat  release,  flame  length,  scorch  height,  and  fireline 
intensity.  The  reaction  intensity,  per  se,  is  not  derived 
from  compensating  factors  and  must  represent  as 
accurately  as  possible  what  it  is  purported  to  be.  So,  the 
experimental  method  of  determining  reaction  zone  size 
and  hence  reaction  intensity  must  be  carefully  defined.  In 
specific  applications,  the  experimental  method  of  deter- 
mining "reaction  intensity"  may  be  significant  for  pre- 
dicting specific  fire  phenomena. 

FIRE  MODEL  GEOMETRIC  EQUATIONS 

Analysis  of  the  geometric  variables,  o,  ft,  and  6,  are 
confined  to  tests  of  the  original  fire  equations  against 
the  expanded  data  base.  The  new  data  are  restricted  to 
the  case  of  no  wind,  zero  (level)  slope,  and  a  minimum  of 
variation  in  fuel  chemistry. 

Rothermel's  equations  are  reduced  to  a  factor-product 
form  for  testing  against  the  new  data  by  the  following 
substitutions  in  equations  27,  36,  37,  and  38  of  table  1.^ 


=  ln(. 


Pp'^1s'"max(''  -  St) 


100 


) 


+  A  (1  -  ln(/?„p) ) 


C,  = 


A  +  1 


a  =  -A/A 


op 


A  fourth  parameter,  C^,  is  introduced  to  test  the  linearity 
of  reaction  intensity  with  loading  and  note  that  in  the 
original  formulation      =  1. 
The  reaction  intensity  equation  becomes 


Ir 


.  /  dW  I      ^C,  „C,  ^CjR  ,C4 


(5) 


where  the  C|  coefficients  are  dependent  on  fuel  size,  o, 
alone.  Anticipating  a  linear  regression  technique,  the  log 
transform  of  equation  5  is: 


In(Ip)  =C^  +  C2ln(/?)  +  C3/J  +  C^\n{6)  +  Hrj^^). 


(6) 


Note  that  in  the  ovendry  limit      =    (a'so  that  in 
Rothermel's  equations      =  1).  Evaluation  of  the  co- 
efficients by  least  square  error  fit  and  comparison  with 
the  original  values  will  test  the  fire  model  formulation. 


Similarly,  the  reaction  velocity,  r',  becomes 


r'  = 


(100  e^^  '\ 
"p  %  J 


(7) 


The  optimum  packing  ratio  ft  is  defined  as  the  packing 
ratio  at  the  maximum  value  of^ the  reaction  velocity. 


Setting 


dV 


dft 


=  0 


we  have  ft^^  =  (1 


op 

whereby  the  optimum  packing  ratio  is  determined  by  the 
regression  coefficients     and  C3. 

In  similar  fashion,  the  propagating  flux  ratio,  |,  is 
expressed  as: 

I       =        ip/lp         =        R(Pp£Qig)  /(-h^). 

Again  substituting  coefficients  in  equation  42  of 
table  1: 

d^  =  [0.0792  +  0.37597  cP-^  -  ln(192  +  7.9095o)] 
dj  =  [0.792  -I-  3.7597cp5] 
the  propagating  flux  ratio  becomes 
I  =  exp  (d^  -1-  d^ft) 


(8) 


or 


ln(Ip/Ip)  =  d,  +  62ft, 


(9) 


where  d.,  and  dj  are  dependent  on  fuel  particle  size,  a. 

Again,  the  coefficients,  d|,  can  be  determined  from  the 
data  and  compared  to  Rothermel's  original  values. 

In  addition  to  the  above,  we  introduce  a  dimensionless 
geometric  parameter— the  product,  oftd,  which  is  the  total 
surface  area  of  all  fuel  particles  per  unit  area  of  fuel  bed. 
We  note  that  heat  transport  out  of  the  combustion  zone, 
heat  absorption  by  fuel  particles,  moisture  mass 
transport  into  and  out  of  fuels,  production  of  volatile 
combustion  products  by  pyrolysis,  etc.,  are  all  related  to 
fuel  surface  area  (further  theoretical  consideration  is 
beyond  the  scope  of  this  paper).  This  product  (oft6)  is 
used  below  in  exploration  of  the  limits  of  combustion. 

The  product  (oft)  has  the  physical  character  of  a  radia- 
tion extinction  coefficient  and  is  introduced  below  in 
speculation  concerning  the  propagating  flux  ratio. 

MOISTURE  DAMPING 

The  empirical  form  of  the  moisture  damping  coeffi- 
cient, r]^,  used  by  Rothermel  is 


(10) 


This  arbitrary  power  series  was  easily  fitted  to  the  limited 
data  available  and  met  the  requirements  that  in  the  dry 
fuel  limit 


^  Note:  the  factor  (in  the  equation  for  ),  1/100,  is  due  to  measurement 
of  fuel  bed  depth  in  centimeters  rather  than  meters. 
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Table  1.— Basic  fire  spread  equations 


Equations 


Units  and  dimension 


Equations 
(Rothermel  1972) 


R  = 


where: 


Ir«1  +      +  ♦=) 


f  = 


'".axWop)*exp[A(1-W,p)l 
(0.0591  +  2.926o-'^"' 

-0.8189 


^     =  0.204O 
A  =  (6.723o°'  -7.27) 
=  1-2.59 

=  0.174S: 


<  -I  .£.1  ) 

V^y     v^y  v^y 


i  = 

=  q0.305U) 


(192  +  7.9095o)~'exp[(0.792  +  3.760o°^)((3  +  0.1)] 

-E 


GO 


C  =  7.47  exp(- 0.87110055) 

B  =  0.1 599  0°^ 

E  =  0.715exp(- 0.01094o) 

Wn  =  WoO-Sr) 

+3  =  5.275/3 -o^dan^^ 

Pb  =  *o'<' 

t  =  exp(  -  4.528/0) 

q  =  581  +  2594M, 


Rate  of  spread,  m/min 
Reaction  intensity,  (kJ/min)/m^ 

Optimum  reaction  velocity,  min'^ 
Maximum  reaction  velocity,  min"' 
Optimum  pact^ing  ratio 

Moisture-damping  coefficient 

Mineral-damping  coefficient 
Propagating  flux  ratio 
Wind  coefficient 


P  = 


Net  fuel  loading,  kg  /  m^ 
Slope  factor 

Ovendry  bulk  density,  kg/nr' 
Effective  heating  number 
Heat  of  preignition,  kJ/kg 
Packing  ratio 


(52) 
(27) 

(38) 
(36) 
(37) 
(39) 

(29) 

(30) 
(42) 
(47) 

(48) 
(49) 
(50) 
(24) 
(51) 
(40) 
(14) 
(12) 
(31) 


Input  parameters  for  basic  equations 


ovendry  fuel  loading,  kg/m? 
fuel  depth,  m 

fuel  particle  surface-area-to-volume  ratio,  cm"^ 
fuel  particle  low  heat  content,  kJ/kg 
ovendry  particle  density,  kg  /  rn' 


fuel  particle  moisture  content. 


S^,  fuel  particle  total  mineral  content. 


gm  moisture 
gm  ovendry  wood 

gm  minerals 


S  ,  fuel  particle  effective  mineral  content, 


gm  ovendry  wood 

gm  silica-free  minerals 


gm  ovendry  wood 
U,  wind  velocity  at  midflame  height,  m/min 
tan  f  slope,  vertical  rise/horizontal  distance 
Mj,  moisture  content  of  extinction,  dimensionless  fraction 
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r,^  =  1  at  M,  =  0 
and  at  the  extinction  limit  of  combustion 

Hm  =  0  at  M,  =  M^. 

The  requirement  that  the  moisture  damping  function  also 
provides  for  extinguishment  at  the  marginal  limits  of 
burning  has  not  worked  satisfactorily  in  application, 
primarily  because  we  have  been  unable  to  find  a  consis- 
tent rationale  for  choosing  a  value  of  the  moisture  of 
extinction,  M^. 

An  alternative  approach  suggested  here  (and 
supported  by  analysis  below)  is  to  separate  extinction 
from  moisture  damping.  A  separate  probability  function 
is  proposed  to  predict  the  marginal  burning  state  of  the 
fire  and  then,  if  it  will  burn,  a  less  complicated  form  of 
the  damping  coefficient  is  needed  to  characterize  fuel 
moisture  effects. 

From  among  several  alternatives,  a  good  fit  to  the  new 
data  is  given  by 

=  expCCgMf)  (11) 

which  meets  the  only  requirement  that  r^^  =  ^  aX  =  0 
for  interpretation  of  the  dry  fuel  limit  in  equation  5  above. 
Then  regression  equation  6  becomes 

ln(Ip)=C^  +  C2ln(/3)  +  C^p  +  C^\n(6)  +  CgM,.  (12) 

The  marginal  burning  probability  function  is  discussed 
in  the  results  section  below.  By  this  rationale  the  mois- 
ture damping  function  will  characterize  the  effects  of  fuel 
moisture  on  Ip  of  burning  fires;  an  independent  function 
will  predict  the  probability  that  it  will  burn  or  not  burn. 

EXPERIMENTAL  SETUP 

Three  size  classes  of  fuel  were  studied:  Populus  spp. 
(poplar)  excelsior  (o  =  81.3  cm-'');  ponderosa  pine,  1/4- 
inch-square  sticks  (o  =  6.30cm-'');  and  ponderosa  pine, 
1/2-inch  sticks  (o  =  3.5  cm-''). 

The  fuels  used  in  these  experiments  were  chosen  to 
minimize  the  effects  of  physiochemical  variations  (fuels 
chosen  were  high-cellulose,  low-ash  woods  and  were 
milled  to  precise  dimensions).  The  physiochemical  char- 
acteristics were  determined  by  standard  techniques  and 
are  tabulated  in  table  2. 

A  series  of  preliminary  fires  were  burned  to  establish 
the  range  of  fire  phenomena  over  the  experimental  range 
of  interest  of  fuel  moisture  and  loading.  These  fires  were 
also  used  to  check  out  instrumentation  and  measure- 
ment techniques.  Preliminary  test  results  are  given  in 
appendix  table  6. 

Table  2.— Physiochemical  characteristics  of  fuels 


Surface/  Particle 
volume  density 
 Fuel  (o),  cm^  ( Pp),  kg/m^ 

Excelsior  (poplar)                      81.3  466 

1/4-Inch  sticks                           6.30  442 
(ponderosa  pine) 

1/2-inch  sticks                           3.15  444 
(ponderosa  pine) 


Table  3  shows  the  matrix  of  variables  examined.  For 
scheduling  and  accounting  purposes,  a  group  number 
was  assigned  to  each  specific  combination  of  o,  p,  and  6 
such  that  for  each  fuel  size,  data  would  be  obtained  for  a 
series  of  groups 

•  at  one  constant  fuel  depth  for  all  packing  ratios; 

•  at  one  constant  packing  ratio  for  all  depths; 

•  at  least  one  constant  load  for  all  depths;  and, 

•  within  practical  limits,  loadings  and  packing  would 

be  duplicated  for  the  different  fuel  sizes. 
Within  each  group  (o,  p,  6),  fires  were  burned  at  various 
fuel  moisture  contents  from  near  ovendry  up  to 
M,  =  0.30.  The  moisture  saturation  limit  of  wood  fiber  is 
approximately  30  percent  of  the  ovendry  wood  weight 
and  is  the  upper  limit  of  woody  fuels  ability  to  absorb 
moisture  from  the  atmosphere  without  immersion  in 
liquid  water. 

Fuels  were  prepared  for  a  series  of  burns  at  a  given 
fuel  moisture  content  by  conditioning  at  controlled  tem- 
perature and  relative  humidity  for  3  to  10  days,  depend- 
ing on  fuel  size  and  target  fuel  moisture.  Prior  to 
conditioning  of  the  excelsior,  the  dry  weight  load  cor- 
rected for  storage  equilibrium  moisture  content  (EMC)  for 
each  burn  was  measured  out  into  baskets.  To  construct 
the  fuel  bed,  excelsior  from  the  conditioning  basket  was 
spread  uniformly  over  the  bed,  fluffed  to  2  or  3  times  the 
planned  fuel  depth  to  remove  tightly  packed  concentra- 
tions, and  lightly  and  uniformly  pressed  down  to  the  final 
fuel  depth  (fig.  1A,  B).  Extraneous  vertical  strands  above 
the  bed  surface  were  then  clipped.  Packing  ratio  for 
excelsior  was  calculated  from  p  =  p6. 

The  first  10  burns  were  loaded  heavily  (assuming  p  - 
32  Ib/ft^)  before  density  measurements  were  completed 
from  fuel  samples. 

For  the  larger  pine  fuels  (1/4-  and  1/2-inch-square 
sticks),  the  packing  ratio  was  controlled  directly  by  con- 
structing each  horizontal  layer  (4  or  2  layers  per  inch  of 
depth,  respectively)  with  the  planned  (solid  volume/total 
volume)  ratio  (figs.  2,  3). 

Raw  data  were  collected  by  an  automatic  digital 
system  (fig.  4)  consisting  of  a  programmable  controller, 
clock,  80-channel  scanner,  digital  multimeter,  thermo- 
couple reference  junctions,  digital  weighing  system, 
35  mm  camera,  and  associated  interfaces  and  sensors. 
(Manual  data  acquisition  consisted  of  visual  measure- 
ments of  the  flame  depth  and  flame  height,  a  photograph 
of  the  fuel  bed  prior  to  ignition,  and  a  tailing  edge  view  of 
the  fire.) 


Effective 

Heat  of  Total  mineral  mineral 

combustion  content  (S^),  content  (Sg), 

(h),  kJ/kg  fraction  fraction 

17  550  0.00347  0.00322 

19  660  .00238  .00224 

19  163  .00367  .00322 
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Table  3. — Fuel  loading  matrix.  First  number  is  the  I.D.  index  to  identify  a  group  of  burns  (series  of  fuel  moistures)  with  constant  loading  and 
geometry.  Second  number  is  the  (ovendry)  fuel  load  in  kglrr^ 


Fuel  bed  Surface  Packing  ratio,  p 
depth  volume   


(cm)  o(cm-^)  0.005  001  002  O04  O08  016  032 


1.27 

81.3 
6.3 
3.15 

3, 

0.050 

9, 

0.201 

2.54 

81.3 
R 

o.o 

3.15 

1, 

0.050 

4, 

.101 

8, 

0.201 

10, 
16, 

.402 
.448 

11, 
20, 
27, 

0.804 
.895 
.895 

22 

1 .790 

33, 

3.580 

5.08 

81.3 
6.3 
3.15 

5, 

.201 

13, 

.448 

17, 

.895 

21, 
28, 

1.790 
1.790 

23, 
31, 

3.580 
3.580 

10.16 

81.3 
6.3 
3.15 

2, 

.201 

6, 

* 

.402 

14, 
24, 

.895 
.895 

18, 
25, 

1.790 
1.790 

29, 

3.580 

32, 

7.161 

34. 

14.322 

20.32 

81.3 
6.3 
3.15 

7. 
12. 

.804 
.895 

15, 

1.790 

19, 
26, 

3.580 
3.580 

30, 

7.161 

'Indicated  individual  special  purpose  fires:  to  test  the  low  moisture  limit  versus  pacl^ing,  etc.  (labeled  Group  #0  in  appendix  table  5). 
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Figure  1.— Excelsior  fuel  bed  for  burn  No.  174(A)  and 
trailing  edge  view  of  the  fire  (B).  o  =  81.3  cnn-'',  p  =  0.04, 
6  =  1.27  cm,  Wq  =  0.201  kg/m2,      =  0.103, 
R  =  0.269  m/nnin,  and  MBI  =  0.6. 
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Figure  2.— Preburn  1/2-inch-square-stick  fuel  bed  (A)  and 
trailing  edge  view  of  fire  (B).  Burn  No.  190.  o  =  3.15 
cm-^  p  =  0.16.  6  =  10.16  cm,      =  7.16  kg/m2_ 
Mj  =  0.165,  R  =  0.50  m/min,  and  MBI  =  1.0. 
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Figure  3.— Preburn  1/4-inch-square-stick  fuel  bed  (A)  and 
trailing  edge  view  of  fire  (B).  Burn  No.  198.  o  =  6.30 
cm-\l3  =  0.08,  6  =  5.08  cnn,       =  1.790  kg/m2, 
=  0.181,  R  =  0.039  m/min,  and  MB!  =  0.9. 
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Figure  4.— Schematic  of  experimental  instrumentation. 
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Immediately  prior  to  ignition,  fuel  bed  samples  were 
bottled  and  sealed  for  later  determination  of  fuel 
moisture  by  xylene  distillation.  Temperature  and  relative 
humidity  in  the  combustion  chamber  were  controlled  at 
approximate  conditions  for  the  desired  EMC  during 
construction  of  the  fuel  bed  and  for  the  duration  of  the 
test  fire. 

Spread  Rate 

The  measurement  of  spread  rate  was  straightforward. 
A  9  row  by  4  column  (fig.  4)  array  of  type  K  (Cr/AI) 
thermocouples  was  suspended  2  to  3  cm  above  the  fuel 
bed  surface  on  a  spacing  of  17.8  cm  (7  in)  between 
columns  and  30.48  cm  (1  ft)  between  rows.  The  TC's 
voltage  was  monitored  at  1.5  to  2  s  cyclic  rate.  When  the 
thermal  EMF  exceeded  10  mv  threshold  (set  to  5  mv  for 
low-intensity  fires)  indicating  arrival  of  the  flame  front, 
the  clock  time,  T  (i,j),  was  recorded  for  each  thermocouple 
(identified  by  the  jth  row  and  ith  column).  The  spread 
rate,  R,  was  calculated  by  least  squares  fit  to  the  equa- 
tion 

where 

T  is  elapsed  time  and  X  is  distance  burned  (fig.  5). 

Reaction  Intensity 

Reaction  intensity  is  defined  as  Ip  =  -h  (dw/dt).  The 
heat  of  combustion,  h,  was  measured  for  each  lot  of  fuel 
by  oxygen  bomb  determination  (table  2).  The  mass  loss, 
dm/dt  (fig.  6),  was  calculated  from  weight  and  time 
measurements  recorded  as  the  fires  progressed  over  the 
last  three  feet  of  the  fuel  bed  (fig.  4).  The  unit  area  weight 
loss  rate,  dw/dt,  is  the  raw  mass  loss  rate  divided  by  the 
combustion  zone  area.  The  area  of  the  combustion  zone 
is  determined  by  the  width  of  the  fuel  bed  and  the  flame 
depth  discussed  below.  The  weight  loss  rates  (listed  in 
column  13  of  appendix  table  5,  are  further  corrected  for 
fuel  moisture  by  1/(1  +  M^)  where     is  the  fractional 
fuel  moisture  from  the  xylene  distillation  determinations. 

Flame  Depth/ Residence  Time 

Rothermel  (1972)  defines  reaction  time,  tp,  as  the 
"time  taken  for  the  fire  front  to  travel  a  distance  equiva- 
lent to  the  depth  of  one  reaction  zone."  Reaction  zone, 
here  called  flame  depth  (D^),  has  never  been  precisely 
defined,  much  less  the  definition  standardized. 

It  is  left  to  the  experimentalist  to  define  this  reaction/ 
combustion/flame  zone  depth  by  his  specific 
measurement  technique.^  Further  difficulty  is  encoun- 
tered in  measuring  the  area  of  the  combustion  zone  (by 
whatever  methods  including  the  present),  because  of  its 
large  inherent  variance  compared  to  other  parameters. 


'The  author  fully  appreciates  the  following  view  expressed  by  F.  A. 
Albini  (private  communication  1981):  "The  problem  here  is  that  the 
idealization  of  the  spreading  fire  as  exhibiting  two  distinct  event 
boundaries,  one  marking  the  onset  of  flaming  combustion  and  one  its 
termination,  is  not  a  wholly  accurate  one.  We  seek  to  preserve  this 
idealization  of  the  process  for  its  power  as  a  modeling  concept,  leaving 
the  experimentalist  the  ill-defined  task  of  fixing  an  operational  definition 
of  its  measurement." 


Anderson  (1964)  and  Rothermel  and  Anderson  (1966) 
calculated  the  horizontal  combustion  zone  depth,  D^,  by 
measuring  the  residence  time,  t,  that  accrued  while 
embedded  thermocouples  remained  above  a  threshold 
temperature  as  the  fire  passed  at  a  speed,  R.  Then 
Dj  =  tR.  Rothermel  (1972)  measured  a  similar  residence 
time  from  the  weight  loss  curves  as  the  combustion  zone 
passed  over  the  leading  edge  of  the  weighing  platform. 

In  the  tests  reported  here,  the  combustion  zone  depth 
is  the  horizontal  distance  from  the  leading  to  trailing 
edge  of  the  base  of  the  flame.  Flame  depth  was  meas- 
ured by  two  direct  methods;  occasionally  these  results 
were  checked  with  results  using  Rothermel's  method 
when  data  and  burning  conditions  were  deemed  appro- 
priate (fig.  7A)— i.e.,  when  no  observable  burnout  existed 
and  the  burning  fuel  did  not  collapse  onto  nor  bridge  the 
weighing  platform  during  transition  of  the  combustion 
zone. 

The  primary  experimental  measurement  of  flame  depth 
(listed  in  column  15  of  appendix  tables  5,  6)  is  the  direct 
measurement  (meter  stick)  of  the  distance  from  leading 
to  trailing  edge  of  the  flames.  Eight  such  measurements 
were  taken  manually  and  averaged  for  each  fire  (inset  fig. 
7B). 

In  addition,  several  flame  probes  were  constructed  of 
linear  arrays  of  thermocouples  for  immersion  in  the  flame 
1  to  2  cm  above  the  fuel  bed  (fig.  4).  The  thermocouple 
temperature  profiles  were  measured  coincident  with  the 
eight  manual  measurements  of  flame  depth  (fig.  7B).  In 
the  example  shown  (burn  number  141),  flame  depth,  Df,  is 
calculated: 

Df  =  (I^i^iV^^ax- 

Flame  depths  measured  in  this  manner  are  consis- 
tently 15  to  25  percent  larger  than  the  manual  measure- 
ments and  invariably  smaller  than  those  calculated  by 
the  weight  loss  technique. 

The  thermocouple  probes  were  incapable  of  resolving 
flame  depths,  smaller  than  5  to  6  cm  (less  than  two  inter- 
vals between  thermocouples)  or  of  measuring  large  zones 
when  the  flame  depth  approached  the  maximum  length 
of  the  probe. 

We  subjectively  observed,  during  the  burning  and 
analysis  of  these  fires,  that  our  experimental  methods 
(Rothermel  1972;  Anderson  1964;  and  present  studies) 
may  overestimate  the  combustion  zone  area  that  con- 
tributes to  propagation  of  the  leading  edge  of  the  fire. 
This  is  particularly  likely  to  happen  in  heavy  fuel  loadings 
where  the  trailing  edge  burnout  is  not  in  a  position  to 
contribute  heat  flux  to  driving  the  fire. 

These  efforts  to  measure  the  total  heat  release  of  the 
combustion  zone  may  provide  (in  most  cases)  adequate 
empirical  relations  for  estimating  flame  lengths  and 
other  fire  effects.  But  the  empirical  calculation  of  a  gross 
"propagating  flux  ratio"  masks  the  physical  process  by 
which  the  propagating  heat  flux  drives  the  fire.  It  should 
be  possible  to  define  or  experimentally  determine  a 
characteristic  path  length  through  (or  flame  depth  into) 
the  combustion  zone  beyond  which  little  or  no  heat  flux 
can  reach  the  unburned  fuel  and  by  which  fire  propa- 
gation could  be  physically  described,  namely,  by  a 
balance  of  driving  heat  flux  and  preignition  heating. 
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Least  So.  Fit 

T  =  18.32851852  +   (26.  5175632  irin/m)  X 

R  square  =  0.999   ;         Spread  Rate  =  0.038  meters/rrinute 
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Figure  5.— Plotted  data  for  rate  of  spread  for  burn  No.  141. 
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Figure  7.— Combustion  zone  depth  from  plotted  weight  loss  rate  (A) 
and  from  flame  probe  temperature  profile  (B)  for  burn  No.  141. 
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This  concept  of  an  effective  combustion  zone  depth 
might  also  be  extended  to  eventually  include  the 
effects  of  wind  on  the  fire. 

Flame  Height 

Flame  heights  are  not  required  in  this  analysis,  but  are 
recorded  here  to  facilitate  associated  fire  modeling 
studies.  Visual/manual  measurements  of  flame  height 
(Hj)  column  14.  appendix  tables  5  and  6,  were  made  simul 
taneously  with  the  flame  depth  measurements  (8  obser- 
vations for  each  fire).  These  are  the  observer's  estimate 
of  the  spacial  and  temporal  average  flame  height  across 
the  flame  zone  for  several  (5  to  10)  seconds.  The  flame 
height  measurement  is  particularly  prone  to  observer 
bias.  The  individual  "official"  experimental  observer 
making  the     and     measurements  was  changed  at 
burn  number  61.  In  addition,  periodically  one  or  two  (of 
six  individuals)  were  asked  to  make  simultaneous  inde- 
pendent Hf  and  D,  measurements.  A  linear  bias  factor 
was  calculated  by  least  square  regression  routine  with  a 
dummy  variable  for  each  observer  for  each  of  the  eight 
observers;  the  resulting  H;  (and  D^)  measurements  listed 
in  appendix  tables  5  and  6  are  normalized  to  the  (D^) 
estimate  of  the  average  observer. 

Another  precautionary  note  concerns  data  taken  of 
flames  by  photography— particularly  by  cameras  with 
automatic  exposure  adjustments.  It  is  obvious,  but  some- 
times forgotten,  that  the  dimensions  of  flames  in 
photographs  are  very  much  dependent  on  the  camera 
exposure  settings.  Longer  exposures  show  higher  peak 
flame  tips  and  wider  flame  bases  than  do  shorter  ex- 
posures. Flame  photographs  were  taken  coincidently 
with  the  eight  measurements  of  each  fire,  but  were  found 
to  be  unusable  for  numerical  data  on  flame  size.  Never- 
theless, the  photos  proved  valuable  for  estimating  the 
general  state  of  combustion. 

Index  of  Marginal  Burning 

It  became  evident  during  the  study  that  the  effect  of 
fuel  moisture  on  burning  rate  followed  a  regular  pattern, 
dependent  primarily  on  fuel  size  class  and  moisture 
content.  The  trend  was  monotone  and  regular  except  that 
the  moisture  content  at  which  the  fires  became  erratic  or 
went  out  varied  widely  and  unpredictably.  Thus  an  ar- 
bitrary "Marginal  Burning  Index"  (MBI)  was  constructed 
with  values  between  -  1  and  1  to  indicate  each  fire's 
proximity  to  extinction.  Fires  that  went  out  (self- 
extinguished)  were  given  an  MBI  equal  to  the  negative 
value  of  the  fraction  of  unburned  fuel  bed  length;  a  fire 
that  burned  out  1  ft  into  the  10-ft-long  conditioned  fuel 
bed  was  given  MBI  =  -0.9,  one  that  burned  out  at  4  ft 
had  MBI  =  -  0.6,  and  so  on.  Secondly,  the  self-sustaining 
fires  were  given  an  MBI  value  equal  to  the  fraction  of  bed 
width,  with  flame  based  attached  to  the  combustion  zone 
as  measured  on  the  end  view  photos  (fig.  1).  These  MBI 
values  are  in  column  16  of  appendix  tables  5  and  6. 

DISCUSSION  OF  RESULTS 

The  study  plan  was  designed  so  that  we  could 
partition  the  data  and  examine  each  independent  vari- 
able, holding  all  others  constant.  Because  additional 


data  are  to  be  collected,  only  a  gross  overview  of  pre- 
liminary results  is  presented  here. 

The  statistical  bias  possibly  introduced  by  the  log 
transforms  was  examined  briefly,  found  negligible,  and 
then  ignored;  that  is,  for 

y  -  Ir.  r',|,r,^,,etc. 

under  some  conditions  the  expectation  value, 

<y>  =  exp(<  In  y>-i-  Sf^^l2) 

where  from  the  linear  regression  (eqn.  6.  9,  12) 

<ln  y>  =  IC.X. 
and       is  the  variance  of  In  y. 

In  y 

In  the  following  discussion  the  expectation 
values  <  y  >  (or  predictions,  y  )  are  estimated  by 
neglecting  the  transform  bias,  thus 

<y  >  =  exp  (<  In  y  >), 

Reaction  Intensity,  Reaction  Velocity 

The  self-sustaining  fires  were  partitioned  by  fuel  size 
class,  o,  then  multiple  linear  regressions  performed  on 
the  log  transformed  form  of  equation  6  (equation  12). 

In  (Ir)  =       +  C^Hp)  +  C3/3  +  C4  In  (d)  -I-  CgM,.  (12) 

The  regression  coefficients  are  listed  in  table  4;  for 
comparison,  the  equivalent  coefficients  of  Rothermel's 
(1972)  formulas  are  also  listed.  The  standard  error  asso- 
ciated with  each  coefficient  is  included  in  parentheses, 
and  the  standard  error  and  r^  of  each  regression  is  also 
listed. 

The  apparently  constant  value  of  C,  for  all  three  sizes 
of  fuel  needs  further  confirmation  and  analysis.  From 
Rothermel's  limited  data,  values  for  this  parameter 
increased  with  increasing  fuel  particle  size  (smaller 
surface/volume,  a).  On  cursory  examination  of  equation  5 
(taking  liberty  against  formal  rigor),  as  the  packing  ratio 
nears  its  optimum  value,  (1  -  C2)/C3,  the  reaction  intensity, 
Ip  =  (Wq/  Pp)e^\  The  new  results  imply  that  the  reaction 
intensity  per  unit  mass  of  dry  fuel  at  the  optimum 
packing  may  be  a  constant,"  independent  of  fuel  size  or 
that  a  free-burning  fire  may  be  correctly  modeled  by  a 
burning  rate  (mass  loss)  that  is  proportional  to  the  total 
fuel  surface  area  (per  unit  area  of  bed)— thus  providing  a 
mechanism  for  handling  fuel  beds  of  mixed  fuel  sizes 
(note  also  that  if  the  above  proves  true,  then       may  be 
determined  from  stationary  crib  fires  with  much  better 
experimental  accuracy).  This  sketchy  speculation  will  be 
examined  and  studied  further  with  test  fires  of  fuels  in 
the  intermediate  sizes,  10<o<50  cm-''. 


*This  result  is  similar  to  that  cited  by  Harmanthy  (1972)  where  room 
and  compartment  fires  with  essentially  unrestricted  ventilation  had  nearly 
constant  (0.0062  kg/m^s)  burning  (massloss)  rates  per  unit  area  of  fuel 
surface.  Our  present  equivalent  value  (from     above)  of  approximately 
0.003  i<g/m^s  is  within  the  lower  range  of  values  from  Harmanthy's 
sources.  Thomas  (1960)  found  similar  values  (0.0021  to  0.0062  kg/m^s) 
for  1-cm-square  sticks  in  stationary  crib  fires. 
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The      coefficient  is  very  close  to  unity  for  the  two 
larger  fuel  sizes  and  reinforces  Rothemnel's  assumption 
that  the  reaction  intensity  is  linear  with  fuel  loading.  In 
the  fine  fuel  (excelsior,  a  =  81.3  cnri-^)  the  fractional 
value  of     (0.65)  indicates  a  weaker  dependence  ( of  Ip 
on  loading)  at  these  loadings  and  packing  ratios.  The 
study  was  designed  to  have  similar  loadings  and  packing 
ratios  for  all  fuel  sizes  within  the  physical  limits  of 
packing  the  fine  fuels.  In  many  of  the  fires  in  fine  fuels 
the  trailing  edge  of  the  combustion  zone  was  not  in  posi- 
tion to  contribute  to  driving  the  leading  edge  of  the  fire 
particularly  in  the  deeper,  heavier  fuel  beds— having  the 
effect  of  decreasing  the  dependence  of  reaction  inten- 
sity on  loading. 

We  plan  to  extend  the  fine  fuel  data  base  to  lighter 
loadings  in  the  near  future.  This  extension  will  include 
loadings  near  20  gm  per  square  meter  and  packing  ratios 
of  0.0005,  which  approach  the  light  loading  limits  found 
in  forest  litter  and  grass  fields. 

In  figure  8  A,  B,  and  C,  the  reaction  intensity  curves 
and  data  (corrected  to  M,  =  0)  are  plotted  and  compared 
with  those  of  Rothermel  (1972).  For  ease  of  comparison 
the  weighted  averages  of  data  at  each  packing  ratio  are 
shown. 

In  figure  9  A,  B,  and  C,  the  reaction  velocity  is  simi- 
larly compared.  The  above  remarks  concerning  the  fine 
excelsior  fuels  are  again  apparent;  in  addition,  we  note 
that  Rothermel  forced  the  reaction  velocity  curve  (for  the 
fine  fuels)  through  the  origin  by  heuristic  argument  to 
obtain  a  maximum  r'  and  p^^.  Figure  9  B  and  C,  show 
significant  differences  in  the  old  and  new  maxima  of  the 
reaction  velocity  curves  for  the  larger  fuel  sizes.  These 
are  the  maxima  that  determine  the  optimum  packing 

ratios,  /J„„.  A  review  of  both  the  old  and  new  raw  data 
'  '  op 

indicates  that  Rothermel's  mass  loss  rates  tend  to  level 
off  with  heavier  loading  (larger  /?)  and,  more  importantly, 
his  residence  time  (or  flame  depth)  increased  signifi- 
cantly with  larger  packing  ratios  (see  discussion  of  flame 
depth  measurements  above). 

Lacking  data  for  light  loading  in  the  fine  fuels,  it  is 
premature  to  draw  new  conclusions.  Also,  data  are 
needed  in  the  intermediate  size  fuels  (10  <  o<  50  cm~  ■") 
before  definitive  improvements  are  made  in  the  present 
fire  equations. 

Propagating  Flux  Ratio 

The  regression  coefficients,  dj,  for  the  propagating  flux 
ratio,  I,  (equation  9),  are  included  in  table  4.  The  regres- 
sion curves,  data,  and  Rothermel's  original  curves  are 
compared  in  figure  10  A,  B,  and  C.  Although  the  coeffi- 
cients in  the  empirical  regression  equations  differ 
considerably,  the  original  Rothermel  curves  are  still 
reasonable  approximations— with  the  possible  exception 


of  underestimating  the  propagation  intensity  in  the  fine 
fuels. 

The  data  shown  in  figure  10  suggest  that  a  linear  form 
of  I  (rather  than  the  exponential)  may  be  a  better  approx- 
imation. The  best  least  squares  fit  of  several  alternative 
forms  (each  including  all  experimental  fires)  is 

I  =  0.017  -  0.317  /3  -I-  0.1119  (o/J) 

with  r2  =  0.91,  standard  error  =  0.048;  the  standard  error 
associated  with  the  coefficient  (0.1119)  of  the  product 
{o(i)  is  0.0026. 

The  product  {oft)  has  the  physical  interpretation  of  a 
radiation  extinction  coefficient  for  the  fuel  bed.  Its  sig- 
nificance as  a  parameter  in  the  propagation  intensity  is 
supported  here.  Tests  of  this  relationship  will  be  ex- 
tended by  experimental  data  in  the  intermediate  size 
fuels  (10<  o<50  cm-1). 

Fuel  Moisture 

We  may  interpret  (-I/C5)  as  a  characteristic  fuel 
moisture,  M^,  dependent  on  fuel  particle  size  (surface/ 
volume  ratio): 

ri^  =  exp(C5Mf)  =  exp(-M,/Mo). 

The  characteristic  moisture      is,  then,  the  fuel  moisture 
at  which  the  reaction  intensity  falls  to  1/e  of  its  value  at 
M,  =  0.  It  can  be  argued  heuristically  that  must 
approach  a  constant  value  in  the  fine  fuel  limit  (o  large) 
and  that      should  have  a  small  positive  value  for  very 
large  fuels  (o  small)  and  that  it  must  fit  the  I/C5  values 
in  table  4;  then 

=  1/(2-1-  2.4  exp  (-0.179o)) 

empirically  described  the  relation  of  M^to  fuel  particle 
size,  a,  for  a  expressed  in  cm'^.  Observe  parenthetically 
that  the  factor  0.179  (cm)  in  the  equation  above  is  a 
characteristic  "fuel  particle  surface  depth"  associated 
with  fuel  moisture  and  fire  propagation,  and  that  it  is 
comparable  to  Frandsen's  (1973)  depth  coefficient 
(1/138  ft  =  0.22  cm)  in  his  effective  heating  equation. 

The  function      is  plotted  in  figure  11.  This  curve 
again  demonstrates  the  need  for  confirming  experimental 
data  for  fuels  in  the  range  of  10<  o<50  cm-''  (fuel 
particles  of  1/8-,  1/16-,  and  1/32-inch  thickness)  near  the 
knee  of  the  M^,  curve. 

Figure  12  is  a  plot  of  the  moisture-damping  curve  with 
the  confidence  bands  for  each  of  the  three  fuel  particle 
sizes.  The  bands  are  calculated  from  the  standard  error 
of  the  coefficient  Cg  alone  and  are  given  here  at  the  67 
percent  confidence  level,  which  is  ±  one  standard  error 
(Draper  and  Smith  1966). 

The  curve  fit  (fig.  12)  includes  all  fires  (however  broken 
up,  erratic,  or  marginal)  that  did  not  go  out. 
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Figure  8.— Reaction  intensities  for  excelsior 
(A),  1/4-inch  fuels  (B),  1/2-inch  fuels  (C). 
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Figure  9.— Reaction  velocities  for  excelsior  (A), 
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Table  4. — Regression  coefficients 


0  =  81.3cm"^ 

0  =  6.3  cm'^ 

0  =  3.1 5  cm"^ 

New 

Standard 

Old 

New 

Standard 

Old 

New 

Standard 

Old 

Coefficients 

value 

error 

value^ 

value 

error 

value^ 

value 

error 

value^ 

1 

11.985 

14.95 

11.81 

20.70 

11.542 

23.91 

2. 

.574 

(0.088) 

1.27 

1.104 

(0.076) 

3.075 

1.143 

(0.090) 

4.591 

^3 

-14.85 

(3.32) 

-49.29 

-3.432 

(1.03) 

-45.93 

-3.828 

(.66) 

-45.06 

C, 

4 

.652 

(.035) 

1.000 

1.053 

(.038) 

1.000 

.923 

(.048) 

1.000 

-1.969 

(.17) 

-2.735 

(.250) 

-3.570 

(.335) 

.82 

.95 

.92 

SE 

.25 

.128 

.156 

'  op 

{-) 

.0056 

.030 

.045 

.037 

.080 

-2.633 

-3.258 

-4.032 

-4.465 

-4.183 

-4.633 

34.07 

(2.02) 

34.69 

7.336 

(.819) 

10.23 

3.413 

(.644) 

7.465 

r2 

.75 

.62 

.39 

^  Equivalent  value  calculated  from  Rothermel  (1972). 
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Figure  11.— Fuel  moisture  characteristic  of  fuel  particle  size. 
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Figure  12.— Moisture-damping  coefficient,  with  expected  error  at  67 
percent  confidence  level  (±  one  standard  error  on  the  coefficient, 
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Extinction 

As  indicated  above,  the  product  oftd  is  the  total  fuel 
particle  surface  area  per  unit  area  of  fuel  bed.  One  might 
expect,  as  surface  area  of  fuel  decreases,  that  fires  will 
be  nnore  difficult  to  burn  and,  conversely,  as  the  fuel 
surface  area  becomes  large,  fires  will  burn  more  easily 
and  burn  at  higher  fuel  moisture.  Indeed,  as  the  experi- 
mental burning  progressed,  a  predictive  rule-of-thumb 
was  developed— fires  will  not  burn  if  the  fuel  moisture 
exceeds: 

Mj>(ln  (2  ap6)  )/4, 
or  if 

op6  <  exp  (4M,)/2. 

From  this  subjective  observation  and  iterative  review  of 
the  experimental  data,  a  family  of  curves, 

of)6  =  2(k^2-3)exp  ( -  kMf)  for  0  <  k  <  20 

was  developed  that  partition  the  experimental  fires  at 
relatively  constant  MBI  (see  Marginal  Burning  Index 
above).  Thus,  an  extinction  index,  n^^,  is  defined: 
In  (8  op6) 

~    M,  +  In  ^2  (15) 

as  a  predictive  indicator  of  marginal  burning.  Large 
values  of  n^  indicate  greater  likelihood  of  good  burning; 
small  values,  poorer  burning.  The  experimental  data,  MBI, 
is  compared  with  n^  in  figure  13.  The  index,  n^,  was  cal- 
culated for  each  fire;  values  ranged  from  2.56  to  13.56. 
The  number  of  fires  was  counted  within  successive  incre- 
ments of  n^,  and  classified  (1)  as  a  "no  burn"  with  zero 
flame;  (2)  by  the  fraction  of  flaming  fire  front  (described 
above),  or  (3)  as  a  "steady  state"  fire  with  100  percent 
continuous  flame  front.  The  average  fraction  of  flaming 
front  was  calculated  (for  the  class  2  fires)  within  each 
increment  of  n^.  Thus,  the  left  most  curve  of  figure  13 
represents  a  probability,  P^,  that  fires  may  burn  (not  go 
out). 

=  -0.885  +  0.3077  n^,  (2.875  <  n^  <  6.125).  (16) 

The  central  curve, 

P^  =  -1.077  +  0.3077  n^,  (3.5  <     <  6.75),  (17) 

represents  the  average  fraction  of  fireline  aflame;  and  the 
right-hand  curve, 

Pgg  =  -  1.538  -I-  0.3077  n^,  (5  <  n^  <  8.25),  (18) 

is  the  probability  of  steady  state  burning  (full,  continuous 
flame  front).  Fuels  will  not  burn  if  nx<  2.87;  if  nx>8.25, 
then  fires  will  burn  well  with  a  hot,  continuous  flame 
front. 

Analysis  of  extinction  lacks  the  rigor  that  will  come 
later  along  with  more  definitive  data;  however,  two  points 


are  illustrated;  first,  that  the  conditions  of  marginal 
burning  are  experimentally  managable  and  not  com- 
pletely unpredictable  (ref.  the  M^  rationale  problem); 
second,  the  fuel  moisture-surface  area  relationship  pro- 
vides a  mechanism  whereby  the  marginal  limits  of 
burning  are  conditional  on  fuel  loading  and  particle  size 
as  well  as  fuel  moisture  content. 

SUMMARY 

The  original  assumption  that  the  reaction  intensity  in 
the  combustion  zone  is  linear  with  loading  is  experi- 
mentally supported  within  the  fuel  sizes,  loadings,  and 
packing  ratios  commonly  encountered  in  wildland  fuels; 
further  confirmation  is  needed  in  light  loadings  of  the 
fine  fuels. 

Rothermel's  empirical  relations  for  the  propagating 
flux  ratio  are  judged  adequate  for  current  use,  though 
perhaps  they  underestimate  the  propagating  flux  in  heavy 
loads  of  fine  fuels.  New  evidence  is  found  to  physically 
relate  the  propagating  flux  ratio  to  the  optical  density 
(aft)  of  the  fuel  bed.  Again,  further  development  is 
needed. 

Experimental  difficulties  in  rigorous  determination  of 
reaction  zone  size  remain  unsolved.  Thus,  physical  in- 
terpretation of  reaction  intensities  should  be  made  with 
caution.  Similarly,  the  causes  and  consequences  of  the 
optimum  packing  ratio  should  be  interpreted  with  cau- 
tion; for  example  (in  the  author's  opinion)  the  widely  held 
idea  that  reaction  intensity  is  limited  by  oxygen  supply 
(indraft)  at  heavier  packing  ratios  is  premature  and  has 
been  overemphasized. 

Confirmation  is  also  needed  for  the  hypothesis  formed 
above  that  the  heat  release  rate  (reaction  intensity)  per 
unit  of  mass  loading  is  constant— independent  of  fuel 
size.  Heretofore,  very  little  experimental  evidence  was 
available  concerning  the  relative  weighting  factors  for  the 
several  size  classes  found  in  wildland  fuels.  If  confirmed, 
the  rationale  is  that  in  free-burning  fires,  the  burning 
mass  loss  rate  will  be  proportional  to  the  total  fuel 
surface  area  per  unit  area  of  fuel  bed. 

The  marginal  limit  of  self-sustained  burning  is  a  func- 
tion of  the  fuel  surface  area  as  well  as  fuel  moisture 
content.  Thus,  an  independent  function  is  proposed  (for 
further  investigation)  as  a  predictive  indicator  of  mar- 
ginal burning.  The  form  of  the  moisture-damping  coeffi- 
cient is  much  simplified  by  thus  removing  the  require- 
ment that  it  provide  for  the  fire  extinction  limit.  After 
removal  of  the  extinction  requirement  from  the  fuel 
moisture-damping  coefficient,  it  is  apparent  that  the 
functional  relation  between  moisture  damping  and  fuel 
particle  size  becomes  tractable. 

Additional  experimental  data  are  needed  in  the  inter- 
mediate fuel  sizes  and  in  the  light  loading  limits  of  fine 
fuels  before  definitive  conclusions  can  be  made. 

An  extension  of  the  study  in  the  light  and  intermediate 
fuels  is  underway. 
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APPENDIX  -  TABULATED  EXPERIMENTAL 
DATA 

In  the  following  data  tables  the  columns  in  left-to-right 
order  are:  1,  tape  file  number;  2,  chronologic  order  of  fire; 
3,  group  number;  4,  5,  6,  ambient  conditions  of  test 
chamber;  7,  fuel  moisture;  8,  surface/volume  ratio  of  fuel 
particles;  9,  packing  ratio  (solid/total  volume);  10,  fuel 
bed  depth;  11,  initial  ovendry  loading;  12,  fire  spread  rate; 
13,  dry  mass  loss  rate  (3-ft-wide  bed);  14,  average  flame 
height;  15,  average  flame  depth;  16,  index  of  flame  con- 
tinuity. 
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Experimental  evidence  from  250  test  fires  confirms  the  general 
formulation  of  Rothermel's  1972  mathematical  model  for  predicting  fire 
spread  in  wildland  fuels.  Numerical  coefficients  are  reevaluated  for  the 
empirical  relations  governed  by  fuel  bed  geometry.  Fuel  moisture  damp- 
ing is  functionally  separated  from  extinction  phenomena.  Hypotheses 
are  proposed  that  propagation  flux  is  related  to  the  optical  density  of  the 
fuel  bed;  that  the  burning  mass  loss  rate  of  different  fuel  size  classes  is 
proportional  to  fuel  particle  surface  area;  and  that  the  marginal  limit  of 
self-sustained  burning  is  dependent  on  the  total  fuel  surface  area  in  the 
bed  as  well  as  the  fuel  moisture  content.  Difficulties  of  measurement 
and  physical  interpretation  of  reaction  intensities  are  discussed. 


KEYWORDS:  fire  behavior,  fire  model,  moisture  damping,  extinction, 
reaction  intensity 


The  Intermountain  Station,  headquartered  in  Ogden, 
Utah,  is  one  of  eight  regional  experiment  stations  charged 
with  providing  scientific  knowledge  to  help  resource 
managers  meet  human  needs  and  protect  forest  and  range 
ecosystems. 

The  Intermountain  Station  includes  the  States  of 
Montana,  Idaho,  Utah,  Nevada,  and  western  Wyoming. 
About  231  million  acres,  or  85  percent,  of  the  land  area  in  the 
Station  territory  are  classified  as  forest  and  rangeland.  These 
lands  include  grasslands,  deserts,  shrublands,  alpine  areas, 
and  well-stocked  forests.  They  supply  fiber  for  forest  in- 
dustries; minerals  for  energy  and  industrial  development;  and 
water  for  domestic  and  industrial  consumption.  They  also 
provide  recreation  opportunities  for  miUions  of  visitors  each 
year. 

Field  programs  and  research  work  units  of  the  Station 
are  maintained  in: 

Boise,  Idaho 

Bozeman,  Montana  (in  cooperation  with  Montana 
State  University) 

Logan,  Utah  (in  cooperation  with  Utah  State 
University) 

Missoula,   Montana  (in  cooperation  with  the 
University  of  Montana) 

Moscow,  Idaho  (in  cooperation  with  the  Univer- 
sity of  Idaho) 

Provo,  Utah  (in  cooperation  with  Brigham  Young 
University) 

Reno,  Nevada  (in  cooperation  with  the  University 
of  Nevada) 


